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Abstract. A linear system of differential equations describing a joint motion of 
elastic porous body and fluid occupying porous space is considered. Although the 
problem is linear, it is very hard to tackle due to the fact that its main differ- 
ential equations involve non-smooth oscillatory coefficients, both big and small, 
under the differentiation operators. The rigorous justification, under various con- 
ditions imposed on physical parameters, is fulfilled for homogenization procedures 
as the dimensionless size of the pores tends to zero, while the porous body is ge- 
ometrically periodic. As the results, we derive Biot's equations of poroelasticity, 
equations of viscoelasticity, or decoupled system consisting of non-isotropic Lame's 
equations and Darcy's system of filtration, depending on ratios between physical 
parameters. The proofs are based on Nguetseng's two-scale convergence method 
of homogenization in periodic structures. 

Key words: Biot's equations, Stokes equations, Lame's equations, two-scale con- 
vergence, homogenization of periodic structures, poroelasticity, viscoelasticity. 



Introduction 

In this article a problem of modelling of small perturbations in elastic de- 
formable medium, perforated by a system of channels (pores) filled with liq- 
uid or gas, is considered. Such media are called elastic porous media and they 
are a rather good approximation to real consolidated grounds. In present- 
day literature, the field of study in mechanics corresponding to these media 
is called poromechanics [Ij. The solid component of such a medium has a 
name of skeleton, and the domain, which is filled with a fluid, is named a 
porous space. The exact mathematical model of elastic porous medium con- 
sists of the classical equations of momentum and mass balance, which are 
stated in Euler variables, of the equations determining stress fields in both 
solid and liquid phases, and of an endowing relation determining behavior 
of the interface between liquid and solid components. The latter relation 
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expresses the fact that the interface is a material surface, which amounts 
to the condition that it consists of the same material particles all the time. 
Denoting by p the density of medium, by v the velocity, by P-^ the stress 
tensor in the liquid component, by the stress tensor in the rigid skeleton, 
and by x the characteristic (indicator) function of porous space, we write the 
fundamental differential equations of the nonlinear model in the form 



where d/dt stands for the material derivative with respect to the time vari- 
able. 

Clearly the above stated original model is a model with an unknown (free) 
boundary. The more precise formulation of the nonlinear problem is not in 
focus of our present work. Instead, we aim to study the problem, linearized at 
the rest state. In continuum mechanics the methods of linearization are de- 
veloped rather deeply. The so obtained linear model is a commonly accepted 
and basic one for description of filtration and seismic acoustics in elastic 
porous media (see, for example, [2l[3|,|l]). Further we refer to this model as 
to model A. In this model the characteristic function of the porous space 
X is a known function for t > 0. It is assumed that this function coincides 
with the characteristic function of the porous space x, given at the initial 
moment. Being written in terms of dimensionless variables, the differential 
equations of the model involve frequently oscillating non-smooth coefficients, 
which have structures of linear combinations of the function x- These coef- 
ficients undergo differentiation with respect to x and besides may be very 
big or very small quantities as compared to the main small parameter e. In 
the model under consideration we define e as the characteristic size of pores 
I divided by the characteristic size L of the entire porous body: 



Denoting by w the dimensionless displacement vector of the continuum 
medium, in terms of dimensionless variables we write the differential equa- 
tions of model A as follows: 



p— = div,,{xP^ + (1 



x)r } + pF, 



— + pdiv^.^ = 





div.P + pF, 



UrP 



P = xa^B{x, 



dw 

'dt 



X)a\^{x, w) 



(g + 7r)I, 



p = -apxdivxw, vr = 



a^(l - x)div^.ty. 



dp 



Op ot 
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Here and further we use the notation 

p = XP/ + (1 - x)Ps, D(x, u) = (1/2) (V,n + (V,it)^) . 

From the purely mathematical point of view, the corresponding initial- 
boundary value problem for model A is well-posed in the sense that it has 
a unique solution belonging to a suitable functional space on any finite tem- 
poral interval. However, in view of possible applications, for example, for 
developing numerical codes, this model is ineffective due to its sophistication 
even if a modern supercomputer is available. Therefore a question of finding 
an effective approximate models is vital. Since the model involves the small 
parameter e, the most natural approach to this question is to derive models 
that would describe hmiting regimes arising as e tends to zero. Such an ap- 
proximation significantly simplifies the original problem and at the same time 
preserves all of its main features. But even this approach is too hard to work 
out, and some additional simplifying assumptions are necessary. In terms of 
geometrical properties of the medium, the most appropriate is to simplify the 
problem postulating that the porous structure is periodic. Further by model 

we will call model A supplemented by this periodicity condition. Thus, 
our main goal now is a derivation of all possible homogenized equations in 
the model . 

The first research with the aim of finding limiting regimes in the case 
when the skeleton was assumed to be an absolutely rigid body was carried 
out by E. Sanchez-Palencia and L. Tartar. E. Sanchez- Palencia [3l Sec. 
7.2] formally obtained Darcy's law of filtration using the method of two- 
scale asymptotic expansions, and L. Tartar [3^, Appendix] mathematically 
rigorously justified the homogenization procedure. Using the same method 
of two-scale expansions J. Keller and R. Burridge |2j derived formally the 
system of Biot's equations [5] from model B^ in the case when the parameter 

was of order e^, and the rest of the coefficients were fixed independent of 
e. It is well-known that the various modifications of Biot's model are bases 
of seismic acoustics problems up-to-date. This fact emphasizes importance 
of comprehensive study of model A and model B^ one more time. J. Keller 
and R. Burridge also considered model B^ under assumption that all the 
physical parameters were fixed independent of e, and formally derived as the 
result a system of equations of viscoelasticity. 

Under the same assumptions as in the article the rigorous justification 
of Biot's model was given by G. Nguetseng ^ and later by A. Mikelic, R. 
P. Gilbert, Th. Clopeaut, and J. L. Ferrin in H El [7|. Also A. Mikelic 
et al derived a system of equations of viscoelasticity, when all the physical 
parameters were fixed independent of e. In these works, Nguetseng's two- 
scale convergence method [SI [ID] was the main method of investigation of the 
model B^ . 
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In the present work by means of the same method we investigate all 
possible limiting regimes in the model B^. This method in rather simple form 
discovers the structure of the weak limit of a sequence {2;'^} as £ \ 0, where 

= u^v'^ and sequences {u^} and {v^} converge as e \ merely weakly, 
but at the same time function has the special structure u'^{x) = u{x/e) 
with u{y) being periodic in y. 

Moreover, Nguetseng's method allows to establish asjTiiptotic expansions 
of a solution of model in the form 

w^{x, t) = {wo{x, t) + ewi (a;, t, x/e) + o{e)) , 

where Wq(x^ t) is a solution of the homogenized (limiting) problem, Wi{x, t, y) 
is a solution of some initial-boundary value problem posed on the generic pe- 
riodic cell of the porous space, and exponent /? is defined by dimensionless 
parameters of the model. Distinct asymptotic behavior of these parameters 
and distinct geometry of the porous space lead to different limiting regimes, 
namely, to various forms of Darcy's law for velocity of liquid component and 
of non-isotropic Lame's equations for displacement of rigid component in 
cases of big parameter a\, also, to various forms of Biot's system in cases of 
small parameter a^, and to different forms of equations of viscoelasticity in 
cases when parameters and a a are 0(1) as e \ 0. For example, in the 
case when 

the velocity v'^ = x'^dw'^ /dt of liquid component and the displacement = 
(1 — x'^)t«^ of rigid skeleton possess the following asymptotic: 

v'{x, t) = e^x' (V {x, t, x/e) + o{e)) , 

u%x, t) =£5(1- ;^=) (^u{x, t) + ell {x, t, x/e) + o{e)) . 

At the same time all equations are determined in a unique way by the given 
physical parameters of the original model and by geometry of the porous 
space. For example, in the case of isolated pores (disconnected porous space) 
the unique limiting regime for any combinations of parameters is a regime 
described by the non-isotropic system of Lame's equations. 

On our opinion, the proposed approach, when the limiting transition 
in all coefficients is fulfilled simultaneously, is the most natural one. We 
emphasize that it is not assumed starting from the original model that the 
fiuid component is inviscid, the porous skeleton is absolutely rigid, or either of 
the components is incompressible. These kinds of properties arise in limiting 
models depending on limiting relations, which involve all parameters of the 
problem. 

The articles [1] and [6], as well as the present one, show in favor of such 
a uniform approach, because they exhibit the situations, when various rates 
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of approach of the small parameter to zero yield distinct homogenized 
equations. Moreover, these equations differ from the homogenized equations 
derived as the limit of model under assumption = 0, imposed even 
before homogenization 

Suppose that all dimensionless parameters of the B'^ depend on the small 
parameter e of the model and there exist limits (finite or infinite) 

limQ;^(£:) = fiQ, \imax{£) = Xq, lima^(£) = Tq, 

£\0 £\0 £\0 

We restrict our consideration by the cases when tq < oo and either of the 
following situations has place. 

(J) yUo = 0, < Ao < OO, 

(//) < yUo < oo, Ao = oo, 

(///) < /io, Ao < oo. 

If Tq = oo then, re-normalizing the displacement vector by setting 

w — s> aT-W, 

we reduce the problem to one of the cases (I)-(III). 

In the present paper we show that in the case (J) the homogenized equa- 
tions have various forms of Biot's system of equations of poroelasticity for 
a two-velocity continuum media or non-isotropic Lam'e's system of equa- 
tions of one-velocity continuum media (for example, for the case of dis- 
connected porous space) (theorem [22]) • In the case (//) the homogenized 
equations are different modifications of Darcy's system of equations of filtra- 
tion for the velocity of the liquid component (and as a first approximation 
the solid component behaves yourself as an absolutely rigid body) and as 
a second approximation - non-isotropic Lam'e's system of equations for the 
re-normalized displacements of the solid component or Biot's system of equa- 
tions of poroelasticity for the re-normalized displacements of the liquid and 
solid components(theorem 12. 3| ). Finally, in the case (///) they are the non- 
local viscoelasticity equations or non-isotropic non-local Lam'e's system of 
equations of one- velocity continuum media (theorem 12.41) . 

§1. Models A and B^ 

1.1. Differential equations, boundary and initial conditions. Let 

a domain Q* of the physical space M.^ be the union of a domain Q* occupied 
with the rigid porous ground and a domain corresponding to hollows 
(pores) in the ground. Domain Q*j is called a porous space and is assumed to 
be filled with liquid or gas. Denote by w*{x*,t*) the displacement vector of 
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the continuum medium (of ground or liquid or gas)) at the point x* G Vt* in 
Euler coordinate system at the moment of time t* > 0. Under an assumption 
that the displacement vector w* is small in which amounts to a case when 
deformations are small, dynamics of rigid phase is described by linear Lame's 
equations and dynamics of fluid or gas is described by Stokes equations. At 
the same time we may set that the velocity vector v*j{x* ,t*) in fluid or gas is 
a partial derivative of the displacement vector with respect to time variable, 
i.e., that v*^{x*,t*) = dw*j:{x*,t*)/dt*. This assumption makes perfect sense 
in description of continuous media in domains, where the characteristic size 
of pores I is very small as compared to the diameter L of the domain Q*, i.e., 
/ = eL and e << 1 (see, for example, [21 [3l H] and recall the observation of 
the linearization procedure for the exact nonlinear model in Introduction). 

In terms of the dimensionless variables, not denoted by the asterisk 
below, 

X* = xL, t* = tr, w*f = WfL, w* = WgL, 1 

^ } (1.1) 

F* = Fg, p*j = ppo, p) = pfpo, p* = pspo, J 

the displacement Wf{x,t) and the pressure p{x,t) of fluid and the displace- 
ment Ws{x,t) of rigid skeleton satisfy the system of Stokes equations 

^rPf-^ = div^P^ + pfF, xe^s,t> 0,(1.2) 

(dvo \ f duo \ 

div^^j -I + a^ofx,^! , xenj, t> 0,(1.3) 

p = —apdiVxWf, X e ilf, t > 0(1.4) 

and the system of Lame's equations 

«rPs^^ = div,.P^ + psF, xeQs,t>0, (1.5) 
= ar^{diVxWs) ■I + axB{x,Ws), xens,t>0. (1.6) 

In f ll.ll) -( |L6l) F = F{x,t) is the given vector of distributed mass forces, L 
is the characteristic macroscopic size - the diameter of the domain r is 
characteristic duration of physical processes, po is mean density of air under 
atmosphere pressure, pf and ps are respectively mean dimensionless densities 
of liquid and rigid phases, correlated with mean density of air, g is the value 
of acceleration of gravity, and po is atmosphere pressure. 

Dimensionless constants aj (z = r, z/, . . .) are defined by the formulas 



L u 2p 

CKj- ^7, Oil/ - , CUn - , 

TLgpo rLgpo 
&Pf 1] 2A 

Lg Lgpo Lgpo 
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where n is the viscosity of fluid or gas, u is the bulk viscosity of fluid or gas, 
A and rj are elastic Lame's constants, and c is a speed of sound in fluid. 

For the unknown functions Wf, p, and Wg, the commonly accepted condi- 
tions of continuity of the displacement field and normal tensions are imposed 
on the interface T := dflf fl d^lg between the two phases (see, for example, 
ISIEISI): 

p/ . n = ■ n, xeT, t>0, (1.8) 
= w', xeT, t> 0. (1.9) 

In (11.81) n is the unit normal vector to F. Note that exactly these conditions 
appear as the result of linearization of the exact nonlinear model. 

Finally, system ( ll.2p -( fL6l) . ( I1.8p - (ll.9p is endowed by giving a displace- 
ment field on S* = dQ and at the moment t = and by giving a velocity 
field at the t = 0. Further without loss of generality, with the aim to sim- 
plify the technical outline, we suppose that these boundary conditions are 
homogeneous. 

1.2. Geometry of porous space. In model A a Lipschitz smoothness 
of the interface between porous space and rigid skeleton is the only restric- 
tion on geometry of porous space. In model the porous medium has 
geometrically periodic structure. Its formal description is as follows [H [11] . 

Firstly a geometric structure inside a pattern unit cell Y = (0, 1)^ is 
defined. Let Yg be a 'solid part' of the cell Y. The 'liquid part' Yf is its open 
complement. Set Y^ := Y^ + k, k & , the translation of Yg on an integer- 
valued vector k. Union of such translations along all A;, Eg := Uy^^^jF^'^ is 
the 1-periodic repetition of Yg all over R'^. Let Ef be the open complement of 
Eg in M^. The following assumptions on geometry of Yf and Ef are accepted. 

(i) Yg is an open connected set of strictly positive measure with a Lipschitz 
boundary, and Yf also has strictly positive measure on Y. 

(ii) Ef and Eg are open sets with C°'^-smooth boundaries. The set Ef is 
locally situated on one side of the boundary dEf, and the set is locally 
situated on one side of the boundary dEg and connected. 

Domains and Q'j are intersections of the domain Q with the sets 
and Rj, where the sets and Rj are periodic domains in R'^ with generic 
cells eYg and eYf of the diameter e, respectively. 

Union Yg U Yf is the closed cube Y = {y eR^, < Vi < I, i = 1,2, 3}, 
and the interface F^ = dQindQ'f is the e-periodic repetition of the boundary 
£7 = edYf n dYg all over u. 

Further by x = X"^ we will denote the characteristic function of the porous 
space. 

For simplicity we accept the following constraint on the domain Q and 
the parameter e. 
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Assumption 1.1. Domain Q is cube, Q := (0,1)^, and quantity 1/e is 
integer, so that Q always contains an integer number of elementary cells Yf. 

Under this assumption, we have 

Xix) = x%x) = x{x/s), (1.10) 

where xiv) is the characteristic function of Yy in Y. 

We say that a porous space is disconnected (isolated pores) if 

7 n = 0. 

1.3. Generalized solutions in models A and B"^. Define the dis- 
placement w{x,t) in the whole domain Q by the formula 

w{x,t) = l''f^ll^ ^^JJ^' (1.11) 

and the pressures p{x,t), q{x,t), and n{x,t) by formulas 

p = -apxdiv^w, TT = -ttr/l - x)diVxW, q = p+ (1-12) 

Thus introduced new unknown functions should satisfy the system 

arP-^ =dw,F + pF, (1.13a) 
dw 

P = xa^B>{x, — ) + (1 - x)axB{x, w) - {q + n)I, (1.13b) 

where p = XPf + 0^x)ps- If xi^) = X%x), then p = p" = X^P/ + (1 - X^)Ps- 
Equations fll.l3p are understood in the sense of distributions theory. They 
involve the both equations (11.21) and (II. 5p in the domains VLf and fi^, respec- 
tively, and the boundary conditions (11.81) and (11.91) on the interface F. There 
are various equivalent in the sense of distributions forms of representation of 
equations (I1.13p . In what follows, it is convenient to write them in the form 
of the integral equality 

^ {arpw .^-pF-iP- xa^Bix, w) : B{x, ^) + (1.14) 

((1 - x)(^xB{x, w)-{q + 7t)I) : D(x, il^)}dxdt = 0, 

where ip{x, t) is an arbitrary smooth test function, such that i/? = dil^/dt = 
at the t = T and i/' = on the boundary St of the domain Qt, = 
n X (0,T). 

In (11.141) by y4 : i? we denote the convolution (or, equivalently, the inner 
tensor product) of two second-rank tensors along the both indexes, i.e., A : 
B = tr (B* o A) = E L=i A^^Bj,. 
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Mathematical model (11.121) . (11.13bp . and (I1.14p is endowed by the initial 
and boundary conditions 

w{x, 0) = 0, ^{x, 0) = 0, xen, (1.15) 

w{x,t) = 0, xeS, t > 0. (1.16) 

Here the assumption, that the boundary and initial conditions are homoge- 
neous, is not essential. 

We will call solution of problem (frT2|) . (]1.13b[) . (11141) -f fLTHll a generalized 
solution in model A and, if fl is supplemented with periodic structure, a 
generalized solution in model B^. 

The following assertion states that problem (fLT2|) . (I1.13bp . (04D - (frT6D 
is well-posed. 

Lemma 1.1. Let the interface T between Qg and Qf be piece-wise continu- 
ously differentiable, parameters pf, ps, O-r, Oi^, a^, ax, and ap be strictly 
positive, and assume that F,dF/dt G L'^^Qt)- 

Then, on any finite temporal interval [0,T], problem (I1.12p . ( ll.lSbl) . 
(I1.14p -( ri.l6p has a unique solution {w,p,7t}, and this solution admits the 
bounds 

Ow dw 

max (y^||(l - x)dtv^-g^{t)\\2,n + V^\\{1 - x)V:, — (t)||2,fi 

d'^w _ dw _ d'^w 

+v^l|-^(^)ll2,n + ^/a^\\xd^v^-g^{t)\\2,n) + y/a^\\xdiVx-^\\2,nT + 

a^it;,, C „,^, ,aF,„ , 



dt^ ""'"^ - dt 
where C depends only on T. 

Due to linearity of the problem, justification of Lemma 11.11 reduces to 
verification of bounds ( I1.17P . These appear by means of differentiation of Eqs. 
(I1.13P with respect to t (note that x and p do not depend on t), multiplication 
of the resulting equation by d'^w/dt'^, and integration by parts. Pressures q 
and 71 are estimated directly from Eqs. (11.120 . 

Further the focus of this article is centered solely on model B^, in which 
coefficients of Eqs. (I1.12P and ( ll.lSp depend continuously on the small pa- 
rametere, p = p^(x) = X^P/ + (l~X^)Ps and g*^, vr^} is a corresponding 

generalized solution. We aim to find out the limiting regimes of the model 
as e \ 0. 

52. Formulation of the main results 
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Suppose additionally that there exist limits (finite or infinite) 
\ima^{e) = z/q, \imap{e) = p^, lima„(e) = 779 . 

£\0 e\0 e\0 

lim — ^ = Ui, lim = pi, hm = Ai, 

2 

hm — ■ — = 771, , lim — - = ri2, hm — = p2. 
In what follows we assume that 

Assumption 2.1. Dimensionless parameters in the model -B^ satisfy restric- 
tions 

/Wo, z/Q, Aq ^ < cx); 0<ro + /ii. (2.1) 

All parameters may take all permitted values. For example, if Tq = 
or p~^ = 0, then all terms in final equations containing these parameters 
disappear. 

The following Theorems I2.1H27^ are the main results of the paper. 

Theorem 2.1. Assume that conditions of Lemma \l.l\ hold and that 
{tu^, g^,/)"^, vr"^} is a generalized solution in model . 
The following assertions are valid: 

(I) If 

Aq < CXD, 

then 

max |||if;^(t)| + ^^x"\^ xw'(t)\ + (1 - x")\^ {t)\\\2,n < If, (2.2) 

Iklh.f^T + lb1l2,f^T + Iklb.^T < If, (2.3) 
where Ip = C\\\F\ + \dF /dt\\\2^Qj, and C is a constant independent of e. 

(II) If 

\q = 00, = 00, < Ai < 00, 
then estimates (12.21) . (12. 3p hold for re-normalized displacements 

with re-normalized parameters 

2 2 2 2 2 

^fj, ^fi 
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(Ill) If 

Xq = OO, Hi < oo, 

then for displacements w'^ hold true estimates (12.21) and under condition 

< OO, (2.4) 

for the pressures and in the liquid component hold true estimates (12.31) . 
// instead of restriction (12. 4p hold true conditions 

0<P2; F = V$, ^, \^\eL'{nr), (2.5) 

then 

max^(||(l - xn'^,{axw%t))h,n + \\x' div,{axw%t))\\2,n < (2.6) 

where iP = C\\\F\ + \d^/dt\ + \dF /dt\\\2,nr andC is a constant independent 
ofe. 

These last estimates imply (12.31) . 

Note, that for the last case {Aq = oo, < oo} we can get same estimates 
(12.21) , (12. 3p and (12. 6p , if instead restrictions (12.50 we consider 

Assumption 2.2. 

F = F'{x,t){l-x') 

and sequences {F^} and {dF^ /dt} are uniformly bounded with respect to e 
in L'^i^T). 

Theorem 2.2. Assume that the hypotheses in Theorem \2. 1\ hold, and 

Ao < OO, /io = 0. 

Then functions w'^ admit an extension from j, = ^1 ^ (0, T) into 
VLt such that the sequence {u^} converges strongly in L'^{Qt) one? weakly in 
L'^{{0,T);W^{n)) to the functions u. At the same time, sequences {w'^}, 
{p^}; o^'iT'd {n^} converge weakly in L'^{VLt) to w, p, q, and n, respec- 

tively. 

The following assertions for these limiting functions hold true: 

(I) If III = OO or the porous space is disconnected (a case of isolated pores), 
then w = u and the functions u, p, q, and n satisfy in the domain Qt the 
following initial- boundary value problem: 

^oPl^ = dtv^i^oA^Q : D(x, u) + B'^div^u + Blq - (g + tt) • 1} + pF, (2.7) 
ot^ 
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— TT + Cq : D(x, u) + a^div^u + a\q = 0, (2.8) 
1 1 dp 

— p H vr + div^u = 0, p + t'oP^^^ = Q, (2.9) 

where p = mpf + (1 — m)ps, m = (x) ■ The symmetric strictly positively 
defined constant fourth-rank tensor Aq, matrices Cq, Bq and Bf and constants 
Qq and al are defined below by formulas (15.301) . (15.321) - (15.341) . 

Differential equations (12. 7p -( I2l8l) are endowed with the homogeneous ini- 
tial and boundary conditions 

Ou 

Tou\t=o = ro — \t=o = 0, xeQ, u{x,t) = 0, x e S. (2.10) 
ot 

(II) If pi < oo, then the weak limits u, , p, q, ir of sequences {u'^}, 
{x^iu^}, {p^}, {q'^}, {vr^} satisfy the initial-boundary value problem consisting 
of the balance of momentum equation 

+ ^oPs(l - m)^ + V(g + tt) - pF = (2.11) 

div^{XoAQ : B{x, u) + B^div^u + -Big}, 

where V = dw^ /dt and Aq, Bq and Bf are the same as in (12. 7p . the continuity 
equations (12. Sp . the continuity equation and the state equation 

— p H 71 + div^w-^ = (m — l)diVxU, p + i^oP^^^ = q, (2-12) 

P* Vo (Jt 

and Darcy's law in the form 

du /"* d'^u 
v = m— + Bi{p^,t-T)-{-V^q + PfF-ToPf-^)ix,T)dT (2.13) 

in the case tq > and pi > 0, Darcy's law in the form 

du 

V = m— + B2{pi) ■ (-Vg + PfF), (2.14) 
in the case tq = and pi > 0, and, finally, Darcy's law in the form 

v = B;- — + (ml-Bs)- / {-Wqix,T)+pfF{x,T))dT, (2.15) 

ot Topf Jq 

in the case tq > and pi = 0. 

This problem is endowed with initial and boundary conditions (I2.10p and 
the boundary condition 

v{x,t) - nix) = 0, xeS, t>0, (2.16) 

for the velocity v of the fluid component. 

In (I2.13l) - (l2.16p n{x) is the unit normal vector to S at a point x & S, 
and matrices Bi{pi,t), B2{pi) and {ml — B^) are defined below by formulas 
dMl), fICTp and (15^6]) . 
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Theorem 2.3. Assume that the hypotheses in Theorem \2. 1\ hold, and that 

Ao = oo. 

(I) If < oo and one of conditions ( 12 ■4p or (12. 5p holds true, then se- 
quences {x'^W^}, {p^} and {g^} converge weakly in L'^^Qt) to , p, and q 
respectively. The functions admit an extension from Vt^ x (0,T) into 
VLt such that the sequence {u^} converges strongly in L'^(flT) and weakly in 
L^{{0,Ty,W^{n)) to zero and 

1) if T'o > ^ and fii > 0, then functions v = dw^ /dt, p and q solve in the 
domain Qt the problem (Fi), where 

v = [ B,{fii,t-r)-{-Vq{x,r)+pfF{x,r))dT, (2.17) 
Jo 

P + J^oP:^^-g^ = q, —-^ + div^v = 0; (2.18) 

^) ^/ To = and pi > 0, then functions v, p and q solve in the domain 
Qj- the problem {F2), where v satisfies Darcy's law in the form 

v = B2{pi)-{-Vq + PfF), (2.19) 

and pressures p and q satisfy equations (12.181) .• 
finally, 

3) if tq > and pi = 0, then functions v, p and q solve in the domain 
Qt the problem (-F3), where v satisfies Darcy's law in the form 

v = (ml-Bs)- / {-Vq{x,T)+pfF{x,T))dT, (2.20) 

Topf Jo 

and pressures p and q satisfy equations (12.181) . 

Problems -F1--F3 are endowed with boundary condition (I2.16p . 

(II) If pi < 00 and conditions (12. 5p hold true, then the sequence {axW^} 
converges strongly in L'^{Qt) and weakly in L^((0,T); 1^2^(i7)) to function u 
and the sequence {tt^} converges weakly in L'^{Qt) to the function n. The 
limiting functions u and n satisfy the boundary value problem in the domain 

n 



= div^iA'o : D(x, u) + B'^div^u + B^q - (g + tt) ■ 1} + pF, (2.21) 

— TT + Co : D(x, u) + a^div^u + a^q = 0, (2.22) 

where the function q is referred to as given. It is defined from the correspond- 
ing of Problems -F1--F3 (the choice of the problem depends on tq and pi). The 
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symmetric strictly positively defined constant fourth-rank tensor A^, matrices 
Cq, Bq and Bf and constants and are defined below by formulas ( I5.30I) . 
fl5.32p - fl5.33p . in which we have rjQ = rj2 and Aq = 1. 

This problem is endowed with the homogeneous boundary conditions. 

(Ill) If fJ'i = oo, Pi^,r]^^ < oo andO < Ai < oo, then there exist weak limits 
, p and 71 of the sequences {a^e^'^x'^'^^} > {p^} ^'^^ {^r^} ^^'^^ « strong limit 
u of the sequence {a^e'^u'^} in L'^{Qt) , which satisfy in Qt the. following 
initial boundary-value problem: 

div^{\ikl : ©(x, u) + B'^div^u + B{p - (p + tt) ■ 1} + pF = 0, 

^ = ^ + ^2(l)-(-Vp + p;F), 

11 , )> (2.23) 

— p H vr + diVxW^ = [ni — l)diVxU, 

^ Pi Vi 

— vr + Cq : ©(x, u) + aldiv^u + a^p = 0. 

Here the symmetric strictly positively defined constant fourth-rank tensor Aq, 
matrices Cq, Bq and Bf and constants Oq and are defined below by formulas 
f l5.30p . f l5.32p - ( I5.34p . in which we have rjQ = rji and Aq = Ai. 

This problem is endowed with the homogeneous initial and boundary con- 
ditions. 

(IV) If fii = oo and Ai = oo, then the corresponding problem for dis- 
placements {a^e~'^w^} has been considered in parts (I)-(II) of the present 
theorem. 

Theorem 2.4. Assume that the hypotheses in Theorem \2. 1\ hold, and that 

< -^0 < oo. 

Then weak limits w, p, q and n of sequences {w'^}, {p^}, {n'^} and {q'^} 
satisfy in fir the following initial-boundary value problem: 



Top- 



dt^ 



+ V(g + 7r) - pF 



diVx{A2 : D(x,— ) +A3 



(x, w) + B^diVxW^- 

I [Aiit — t) : D(x, w{x, t)) + B^{t — T)diVxW{x, T))dT). 
Jo 

1 



(2.24) 



-p + mdiVxW 



/ {C2{t — t) : D(x, w{x, r)) + a2{t — T)diVxW{x, T))dT. 
Jo 



(2.25) 
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— 71 + (1 — m)diVxW 



— / (^C3{t — t) : 3{x,w{x,t)) + a3{t — T)diVxW{x,T))dT, 
Jo 



(2.26) 



Here A2, A2 and A4 are fourth-rank tensors, B4, B^, C2 and C3 are 
matrices and 02 and 03 are scalars. The exact expressions for these objects 
are given below by formulas f l7.26p - (17.311) . 

The problem is supplemented by the homogeneous initial and boundary 
conditions, (see in Eqs. (I2.10p in Theorem \2.^) . 

If the porous space is connected then A2 is strictly positively defined sym- 
metric tensor. 

If the porous space is disconnected, which is a case of isolated pores, then 
A2 = and the system (I2.24p degenerates into nonlocal anisotropic Lam'e's 
system with strictly positively defined and symmetric tensor A3. 

§3. Preliminaries 

3.1. Two-scale convergence. Justification of Theorems 12.1112.41 relies 
on systematic use of the method of two-scale convergence, which had been 
proposed by G. Nguetseng [8] and has been applied recently to a wide range 
of homogenization problems (see, for example, the survey pDj). 



Definition 3.1. A sequence {f'^} C L'^{Qt) is said to be two-scale convergent 
to a limit ip G L^{Qt x Y) if and only if for any 1-periodic in y function 
a = a{x, t, y) the limiting relation 

lim / (f^{x,t)a {x,t,x/e) dxdt = / / Lp(x,t,y)a(x,t,y)dydxdt (3.1) 
holds. 

Existence and main properties of weakly convergent sequences are estab- 
lished by the following fundamental theorem [SJ ilOj : 

Theorem 3.1. (Nguetseng's theorem) 

1. Any bounded in L'^{Q) sequence contains a subsequence, two-scale 
convergent to some limit ip G L'^{VLt x Y). 

2. Let sequences {^^} and {e'Wx'-P^} be uniformly bounded in L'^{Qt)- Then 
there exist a 1-periodic in y function ip = ip{x,t,y) and a subsequence {<p'^} 
such that {p, Vyip G L'^{Qt x Y), and ipf^ and eVx'^^ two-scale converge to 
and Vyip, respectively. 

3. Let sequences {^^} and {V^^v^^} be bounded in L'^[Q). Then there exist 
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functions (f G L'^{Qt) ond ip G L'^{Qt x Y) and a subsequence from {v^^} such 
that ip is 1-periodic in y, Vyip G L'^{Qt x Y), and if'^ and Vxf^ two-scale 
converge to ip and Vxf{x,t) + V ytpi^x^tjy) , respectively. 

Corollary 3.1. Let a G LF'iX) and a'^{x) := a{x/e). Assume that a se- 
quence {^^} C L'^{VLt) two-scale converges to (f E L'^{Qt x Y)- Then the 
sequence a'^ip'^ two-scale converges to a(p. 

3.2. An extension lemma. The typical difficulty in homogenization 
problems while passing to a limit in Model 5^ as e \ arises because of the 
fact that the bounds on the gradient of displacement VxW^ may be distinct in 
liquid and rigid phases. The classical approach in overcoming this difficulty 
consists of constructing of extension to the whole Q of the displacement field 
defined merely on Qg. The following lemma is valid due to the well-known 
results from [T^l [13]. We formulate it in appropriate for us form: 

Lemma 3.1. Suppose that assumptions of Sec. 1.2 on geometry of periodic 
structure hold, ip'^ G W2{^1) and = on = dQ^ fl dQ in the trace 
sense. Then there exists a function G 1^2^ (fi) such that its restriction on 
the sub-domain fll coincide with ip'^ , i.e., 

(1 - x'{x)){a'{x) - = 0, xen, (3.2) 

and, moreover, the estimate 

h'hn < C\m\2,ni, W^.^y'hn < CWWx^hMi (3-3) 

hold true, where the constant C depends only on geometry Y and does not 
depend on e. 

3.3. Priedrichs— Poincare's inequality in periodic structure. The 

following lemma was proved by L. Tartar in [31 Appendix]. It specifies 
Friedrichs-Poincare's inequality for e-periodic structure. 

Lemma 3.2. Suppose that assumptions on the geometry of f2j hold true. 

o 

Then for any function ip (fij) the inequality 

I \'p\'^dx < Ce^ I \Vx'^?dx (3.4) 

holds true with some constant C, independent of e. 

3.4. Some notation. Further we denote 
1) 
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{^)n= / ^dx, {ip)n^ = / ipdxdt. 
Jn Jvlt 

2) If a and b are two vectors then the matrix a ® b is defined by the formula 

(a b) • c = a(b ■ c) 

for any vector c. 

3) If B and C are two matrices, then i? ® C is a forth-rank tensor such 
that its convolution with any matrix A is defined by the formula 

{B®C):A = B{C : A) 

. 4) By P-' we denote the 3 x 3-matrix with just one non-vanishing entry, 
which is equal to one and stands in the z-th row and the j-th column. 
5) We also introduce 

X' = + in = ^(e. ® e, + ® e^), 
where (61,62,63) are the standard Cartesian basis vectors. 

§4. Proof of Theorem [211] 

4.1. Let Ao < 00. 

If restriction tq > holds, then (12.21) follow from the lemma 11.11 
Let p^. < 00 and r^o < 00. Then pressures and vr^ are bounded from 
equations (11.121) with the help of (I1.17p . The pressure is bounded from 

=p' - a,.X^div^-^ (4-1) 

with the help of ffrTTj) . 

If j9* = 00, then estimate (12.31) for the sum of pressures (g"^ + vr^) follows 
from the basic integral identity (11.141) and estimates (I1.17P as an estimate 
for the corresponding functional, if we re-normalize the pressures {q^ + vr"^) 
such that 

/ {q'{x,t) +7r'{x,t))dx = 0. 
Jq 

. Indeed, the basic integral identity (I1.14p and estimates (11.171) imply 



Jq 



2,n- 



Choosing now if) such that (g^ + vr^) = div^.i/' we get the desired estimate for 
the sum of pressures (g^ + tt^). Such a choice is always possible (see [I^), if 
we put 

-0 = V(/)+'0o, div^^Vo = 0, Av? = gHvr^ v'lan = 0, (Vv2+V'o)|an = 0. 
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Note that the re-normahzation of the pressures (g^+vr"^) transforms continuity 
and state equations fll.l2p for pressures into 

^f + x'dw^,w' = --(3'x', (4.2) 
ap m 

^-TT^ + (1 - x^)div.tf;^ = (1 - x')ll (4.3) 



q' = f + -% + x'i}. (4.4) 

ap at ■' 



where 



= ((1 - xndiv,w')n, mi) = {q')n, (1 - m)^^ = - — {q')n + P'- 

The case ?7o = oo is considered in the same way. Note that for all sit- 
uations the basic integral identity fll.141) permits to bound only the sum 
(g'^ + TT^). But thanks to the property that the product of these two functions 
is equal to zero, it is enough to get bounds for each of these functions. The 
pressure is bounded from the state equation (14. 4p . if we substitute the 
term {a,j / ap)dp'^ / dt from the state equation (14. 2 p and use estimate (I1.17p . 

Estimation of in the case Tq = is not simple, and we outline it in 
more detail. 

Let /ii > and tq = 0. As usual, we obtain the basic estimates if we 
multiply equations for w'^ by dw^ / dt and then integrate by parts all obtained 
terms. The only one term F ■ dw^ /dt heeds additional consideration here. 
First of all, on the strength of Lemma 13.11 we construct an extension of 
the function from fi^ into fij such that u'^ = in fi^, ■u'^ G PV2^(f2) and 

\\u'\\2,n < C\\V.u'h,n < -^ll(l - x')v^V,io^||2,c. 



After that we estimate ||'i«^||2,n with the help of Friedrichs-Poincare's 
inequality in periodic structure (lemma [3. 2p for the difference {u'^ — w^): 



_ 1 

< \W\\2,n + Ce\\\/^u^\\2,n + C{ea^,'')\\x^ ^/a^/S ^w^\\2,q 

C -i 
< ^11(1 - x")V^x^x'w'\\2,^ + C{ea^'')\\x'^^V^w'\\2,n- 

Next we pass the derivative with respect to time from dw'^ /dt to p'^F and 
bound all positive terms (including the term ayX^diVxdw'^ / dt) in a usual way 
with the help of Holder and Grownwall's inequalities. 
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The rest of the proof is the same as for the case tq > if we use a 
consequence of f l 1.171) : 

max ar||^T^(i)||2,fi < C. 

0<t<T dv 

4.2. The proof of this part of the theorem is obvious, because the re- 
normahzation reduces this case to the case of /ii = 1 and tq = 0, which 
has been already considered. 

4.3. Let Ao = C)0, /ii < oo and conditions (12. 5p hold true. It is obvious that 
estimates (12.21) are still valid. 

The desired estimates (12. 6p follow from the basic equations for a\W^ in 
the same way as in the case of estimates (12. 2p . The main difference here is 
in the term p'^F ■ a\dw^ / dt^ which now transforms to 

T = pfF- ax-g^ + {pf - P/)(l - X')F ■ (^x-gf- 
The integral of first term in T transforms as 

Of / V<P ■ ax— — dxar = —pf / / <Pa\aWx^ — dxdT 
Jo Jn dr Jo Jn dr 

/"* f dw^ dw^ 

= -pf 1 / (X^ ■ $«Adiv^.— — + (1 - x^) ■ ^axdiv^——)dxdT 



'pf 

pf ' 



j ix^ ■ $a;AdiVj,t(;^ + (1 - X^) ■ ^axdi\xu'')dx+ 
Jn 

[ [ ix' ■ ^rttAdiv^iu^ + (1 - x') ■ <l>raxdiyxu')dxdT 
Jo Jn 

and is bounded with the help of terms 

/ {x'{apax^){axdiyxw')^ + (1 - x')\ax^xu'\'^)dx, 
Jn 

which appear in the basic identity after using the continuity equations. 
The integral of the second term in T is bounded with the help of the term 

{{l~xn\(yxVxU'\^)n 

in the same way as before. 

Estimates (12. 3p follow now from (12.60 . Here as before the sum of pressures 
(g^ + 7r^) is bounded from the basic integral identity (I1.14p as a corresponding 
functional, and the pressure is bounded from the equation (14. ip due to 
the bound (11.170 for the divergence of the velocity of the liquid component 
X'^diVx{dw'^ /dt). 

If instead of conditions (12.50 one has condition (12.40 . then bounds (12.30 
for pressures and follow from equations (11.120 and (14. ip and bounds 
fOT]) . Note that in this case (3' = 0. □ 
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§5. Proof of Theorem [Ml 



5.1. Weak and two-scale limits of sequences of displacement 
and pressures. On the strength of Theorem 12.11 the sequences {p^}, {q'^}, 
{tt^} and {w'^} are uniformly in e bounded in L'^{Qt)- Hence there exist a 
subsequence of small parameters {e > 0} and functions p, q, n and w such 
that 

p'^ — > p, — > g, n'^ ^ n, w'^ ^ w (5.1) 

weakly in L'^iVLx) as e \ 0. 

Due to Lemma [3.11 there is a function G L°°(0, T; VF2^(f2)) such that 
u'^ = in fig X (0, T), and the family {u^} is uniformly in e bounded in 
L°°(0, T; W2{fi)). Therefore it is possible to extract a subsequence of {e > 0} 
such that 

weakly in L'^{0, T; W^{n)) (5.2) 

as £ \ 0. Moreover, 

X'cy'^m^') - 0. (5.3) 
Relabelling if necessary, we assume that the sequences converge them- 
selves. 

On the strength of Nguetseng's theorem, there exist 1-periodic in y func- 
tions P{x,t,y), Il{x,t,y), Q{x,t,y), W{x,t,y) and U{x,t,y) such that 
the sequences {p^}, {vr^}, {g^}, {W^} and {V^H^} two-scale converge to 
P{x, t, y), Il{x, t, y), Q{x, t, y), W{x, t, y) and V^u + WyU^x, t, y), respec- 
tively. 

Note that the sequence {div^-io^} weakly converges to div^iu and u G 

o 

L'^{0,T;W2 (fi))- Last assertion for disconnected porous space follows from 

o 

inclusion u'^ G ^^^(0, T; W2 (fl)) and for the connected porous space it follows 
from the Friedrichs-Poincare's inequality for u'^ in the e-layer of the bound- 
ary S and from convergence of sequence to u strongly in L'^{Qt) and 
weakly in L^{{0,Ty,W^{n)). 

5.2. Micro- and macroscopic equations I. 

Lemma 5.1. For all x E Q and y E Y weak and two-scale limits of the 
sequences {p'^}, {tt"^}, {q'^}, {W^}, and {u'^} satisfy the relations 



P = ixp, Q = ixq; 




(5.4) 


+ (1 - x){div^u + diVyU) = 7, 




(5.5) 


diVyW = 0; 




(5.6) 


W = x{y)w + {1 - x)u; 




(5.7) 


q = p + -i^QP~^%+lfm] 




(5.8) 


yP + div^w = (1 — m)diVxU + {divy 




(5.9) 


^vr + (1 - m)div^u + {diVyU)Y, = ( 


1 - m)7„ 


(5.10) 
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where 

(3= f {diVyU)Ysdx, m-ff = (1 - m)-fs = - — {q)n + P, 

Jn Vo 

ifp* + Vo = oo and /3 = 7/ = 7^ = 0, if p^: +r]o < 00. 

Proof. In order to prove Eq. ( 15. 4p . into Eq. ( 11. 14^ insert a test function 
i/)'^ = eil^ (x, t, aj/e), where ^) v) is an arbitrary 1-periodic and finite on 
Yf function in y. Passing to the limit as e \ 0, we get 

WyQ{x,t,y) = Q, yeYf. (5.11) 

The weak and two-scale hmiting passage in Eq. ( 14.41) yield that Eq. ( 15.81) 
and the equation 

7/ BP 

= ^ + + 

hold. Taking into account Eq. ( 15.111) and ( 15.12p we get 

VyP{x,t,y) = Q, yeYf. 
Next, fulfilling the two-scale limiting passage in equalities 

(1 - x')p' = 0, (1 - x')q' = 

we arrive at 

{i-x)P = o, (i-x)Q = o, 

which along with Eqs. ( 15. lip and ( 15.12p justifies Eq. ( 15.40 . 

Eqs. ( 15.51) . (15.61) . (15.91) . and (15.10p appear as the results of two-scale 
limiting passages in Eqs. (14. 2p - (14. 4p with the proper test functions being 
involved. Thus, for example, Eq. (15. 9p arises if we represent Eq. (14. 2p in the 
form 

/ + div.tt;^ = (1 - x')dw.u' - -P'x', (5.13) 

m 

multiply by an arbitrary function, independent of the "fast" variable x/e, 
and then pass to the limit as e \ 0. Eq. ( 15.10p is derived quite similarly. 
In order to prove Eq. ( 15.60 . it is sufficient to consider the two-scale limiting 
relations in Eq. ( 15.130 as e \ with the test functions eip (x/e) h{x,t), 
where ip and h are arbitrary smooth test functions. In order to prove Eq. 
(15.71) it is sufficient to consider the two-scale limiting relations in 

{1 - X'){w' - W") = 0. 

□ 
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Lemma 5.2. For all {x,t) G Qt the relations 

XoAyU = VyU, y e (5.14) 

(XoB(y, [/) - n • I + XoB(x, u) + —q-l)-n = 0, y e-f, (5.15) 

m 

hold true. Here n is a unit normal to 7. 

Proof. Substituting a test function of the form 1/?^ = eip {x,t,x/e), where 
il}{x, t, y) is an arbitrary 1-periodic in y function vanishing on the boundary 
S, into Eq. fll.i4p . and passing to the hmit as £ \ 0, we arrive at the following 
microscopic relation on the cell Y: 

div,{Ao(l - x)my, U) + D(x, u)) - (n + Igx) ■ 1} = 0, (5.16) 

which is clearly equivalent to Eqs. fl5.14p and (15.150 in view of Eqs. (l5.4p . □ 

Lemma 5.3. Let p = mpf + (1 — rn)ps, V = xdw/dt, and v = {V)y- Then 
for all < tq < 00 the quadruple of functions {u,v,q,7T} satisfies in Qt the 
system of macroscopic equations 

dv d'^u 
^oPf-Q^ + rops{l - m)-^ - pF = (5.17) 

div,{Xo{{l - m)B{x, u) + {B{y, C/))yJ - (g + vr) ■ I}. 

Proof. Eqs. f l5.17p arise as the limit of Eqs. fll.l4p with test functions being 
finite in Qj' and independent of e. □ 

5.3. Micro- and macroscopic equations II. 

Lemma 5.4. If pi = 00 then the weak and two-scale limits of {u"^} and 
{w'^} coincide. 

Proof. In order to verify, it is sufficient to consider the difference (it^ — w^) 
and apply Friedrichs-Poincar'e's inequality, just like in the proof of Theorem 
O □ 

Lemma 5.5. Let pi < 00. Then the weak and two-scale limits of {q"^} and 
{x^w^} satisfy the microscopic relations 

dV 

ToPf-g^- PfF = piAyV -VyR-V^q, yeYf, (5.18) 

V = ^, ye^ (5.19) 

in the case pi > 0, and relations 
dV 

roPf— = -VyR-V^q + PfF, yeVf, (5.20) 
22 
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)-n = 0, ye^ (5.21) 

in the case /ii = 0. 

In Eq. (15.211) n is the unit normal to 7. 

Proof. Differential equations (I5.18P and (15.201) follow as e \ from integral 

equality (11.141) with the test function if) = ifi{xe~^) ■ h{x,t), where cp is 

solenoidal and finite in Yf. 

Boundary conditions (15.191) are the consequences of the two-scale conver- 
1 i_ 

gence of {a^VxW'^} to the function filVyWi^x, t, y). On the strength of this 
convergence, the function VyW{x., t, y) is L^-integrable in Y. The boundary 
conditions (15.211) follow from Eq. (15. 6p . □ 

Lemma 5.6. // the porous space is disconnected, which is the case of iso- 
lated pores, then the weak and two-scale limits of sequences {u"^} and {w^} 
coincide. 

Proof. Indeed, in the case < /Ui < 00 the systems of equations (15. 6p . 
(EUD, and (15A911 . or (ESI), (^M), and (ICTl) have the unique solution V = 
du/dt. □ 

5.4. Homogenized equations I. 

Lemma 5.7. // /ii = 00 or the porous space is disconnected then w = u 
and the weak limits u, p, q, and tt satisfy in Qt the initial-boundary value 
problem 

= diVx{\oK ■ K)(a;, u) + B'^diu^u + Blq - (g + tt) ■ 1} + pF, (5.22) 

— TT + Cn : D(x, u) + andiv^u + a\q + aUq)^. = 0, (5.23) 
^0 

dp 1 1 
g = p + z/op~^— + 7/m, — p H n + div^u = (1 - m)'js - P, (5-24) 

where the symmetric strictly positively defined constant fourth-rank tensor 
Aq, matrices Cq, Bq and Bl and constants a^, a\ and are defined below by 
formulas f lCTD . f 02D - f OD . 

Differential equations (15.221) are endowed with homogeneous initial and 
boundary conditions 

du 

tqu{x,0) = tq — {x,0) = 0, xeQ, u{x,t) = 0, X e S, t > 0. 

(5.25) 

Remark 5.1. In what follows we can neglect terms 02(^)0, 7/ and 7^ because 
all pressures are defined up to arbitrary functions of time. 
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Proof. In the first place let us notice that u = w due to Lemmas 15.41 and 



The homogenized equations (15.221) follow from the macroscopic equations 
f l5.17p . after we insert in them the expression 

{B{y, U))y, = Al : 3{x, u) + B'^dw.u + B^q. 

In turn, this expression follows by virtue of solutions of Eqs. (15.51) and (15.161) 
on the pattern cell Yg. Indeed, setting 

3 1 
U=J2 U'Ky)D^, + C/o(?/)div.n + -{U,{y){q - {q)n) + U2{y){q)ci), 



where 



3 1 

n = Ao V ii''{y)D,, + Uo{y)dw^u + -(ni(2/)(g - {q)n) + My){q)i 



(5.26) 

(5.27) 
(5.28) 



we arrive at the following periodic-boundary value problems in K,: 

div,{(i - x)my, u^n + j'n - ■ i} = o, 

^n*. + (1 _ x)divyU'' = 0; 

div,{Ao(l - xMy, Uo) - Ho ■ 1} = 0, 
^no + (l-x)(div,C/o + l)=0; 

div,{Ao(l - xMy, Ui) - (Hi + x) ■ 1} = 0, 
+ (1 - x)divyU,) = 0. 

div,{Ao(i - xMy, U2) - (Hs + x) • 1} = 0, 

1 (1 - y) m ) (5.29) 

-n2 + (l-x)div,C/2-)-^((div,t/2)y. + -) = 0. ^ 
m (1 - m) " r]Q 

Note that for p^, + 770 = 00 

3 

(3 = {{diYyU)n)Ys = ^{diYyU'^)YADij)n + (divj,t/o)n(div,ii)o+ 

— (diVj,C/i)y,((g - {q)^))n + — (divyf/2)y, (g)f7 = — (diVj^C/2)y, 
m mm 

due to homogeneous boundary conditions for u{x,t). 

On the strength of the assumptions on the geometry of the pattern "liq- 
uid" cell Kj, problems (15.261) - (??) have unique solution, up to an arbitrary 
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constant vector. In order to discard the arbitrary constant vectors we de- 
mand 

{U'')y. = (C/o)n = {Ui)y^ = {U,)y^ = 0. 

Thus 



J2 J'' ® J'' + K K\ = Y,{{^-x)D{y,W'))Y®f^. (5.30) 

i,j=l i,j=l 



Symmetry of the tensor Aq follows from symmetry of the tensor And 
symmetry of the latter one follows from the equality 



{Biy,U'^))y^ (5.31) 

-my, U'^) : D(y, U'^))y^ + ^Wm'\ 

Vo 

which appears by means of multiplication of Eq. (I5.26P for U^^ by U^'' and 
by integration by parts. 

This equality also implies positive definiteness of the tensor Aq. Indeed, 
let ( be an arbitrary symmetric matrix. Setting 

3 3 
i,j=l i,j=l 

and taking into account Eq. (15.311) we get 

(D(y,Z))y^ : C = -(E)(y,Z) : B{y,Z))y^ - 

Vo 

This equality and the definition of the tensor Aq give 

(A^ c) : C = my, + C) : my, + C))n + -n'- 

Now the strict positive definiteness of the tensor Aq follows from the equality 
immediately above and the geometry of the elementary cell Yg. Namely, let 
for some function ( such that C • C = 1 be (Aq :():( = 0. But then 
we have (©(y, 1^) + () = 0, which is possible iff Z is a linear function in y. 
On the other hand, all linear periodic functions on Yg are constant. Finally, 
the normalization condition (U^^) = yields that Z = 0. However, this is 
impossible because the functions are linearly independent. 

Finally, Eqs. (I5.23P and (15.241) for the pressures follow from Eqs. (15. 8p - 
(ICTj) and 

= my, C/o))y„ Bl = i-(D(y, t/i))y„ (5.32) 
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3 

Q = J2 {dwyU'')Yj'', = 1 - m + (div,C/o)y,, (5.33) 

1 1 777 

= -(div,t/i)y,, 4 = -(- - (div,t/i)yj. (5.34) 



m m rjQ 



□ 



5.5. Homogenized equations II. 

Let /ii < oo. In the same manner as above, we verify that the hmit u of 
the sequence {W^} satisfies the initial-boundary value problem hkes fl5.22p - 
(I5.24p . The main difference here that, in general, the weak limit w of the 
sequence {W^^} differs from u. More precisely, the following statement is 
true. 

Lemma 5.8. // yUi < oo then the weak limits u, , p, q, and vr of the 
sequences {u'^}, {x^^^}; {p^}; {O'^}; o,'>T'd {vr^} satisfy the initial-boundary 
value problem in Qt, consisting of the balance of momentum equation 

^o(p/^ + Psil + V(g + vr) - pF = (5.35) 

diVx{XoA^Q : D(x, u) + B^div^u + B{q}, 

where V = dw^ /dt and Kq, Bq and Bl are the same as in (12. 7p . the continuity 
equations (15.120 . the equations 

^dp 1 1 , f / 

p + uqp = q, — p H vr + diVxW-' = (m — IjdiVxU, (5.36) 

ot p^ r]o 

and Darcy's law in the form 
du 



V = m—— 

dt 



+ j Bi{i2i,t-T)-{-Vxq + PfF-ToPf-^)ix,T)dT (5.37) 

in the case o/tq > and pi > 0, Darcy's law in the form 

du 

v = m— + B2{pi)-{-Vxq + PfF) (5.38) 
in the case of tq = and fii > 0, and, finally, Darcy's law in the form 

v = B:,- — + (ml-Sg)- / (- V,g(a;, r) + p^F(a;, r))rfr (5.39) 

Ct TqPj Jq 

in the case of tq > and pi = 0. The problem is supplemented by boundary 
and initial conditions (I5.24p for the displacement u of the rigid component 
and by the boundary condition 

v{x, t) ■ n{x) = 0, {x, t) eS = dn, t>0, (5.40) 

for the velocity v of the liquid component. In Eqs. (15.370 - (15. 40p n{x) is the 
unit normal vector to S at a point x E S, and matrices Bi{fj,i,t), B2{pi), 
and i?3 are given below by Eqs. (15.420 - (15. 46p . 
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Proof. The derivation of Eq. fl5.40p is standard [3]. The homogenized 
equations of balance of momentum and balance of mass derive exactly as 
fl5.22p - fl5.23p . For example, to get Eq. fl5.36p we just expressed div^'u; 
in a sum of Eqs. (15.90 and (15.100 using Eq. (15. 7p after homogenization: 
w = + (1 — m)u. Therefore we omit the relevant proofs now and focus 
only on derivation of homogenized equations for the velocity v in the form 
of Darcy's laws. 

a) If /ii > and tq > 0, then the solution of the system of microscopic 
equations (15. 6p . (15.180 . and fl5.19p . provided with the homogeneous initial 
data, is given by formula 

^ = ^ + B({y,t-r)-{-V.q + PfF-ToPf-^){x,T)dT, 
in which 

3 

B{(2/,t) = J]V'(t/,t)®e„ 

i=l 

and the functions V\y, t) are defined by virtue of the periodic initial-boundary 
value problem 

^oP/^ - /iiAV* + Vg^ = 0, dwyV = 0, yeYj,t>oA 
V' = 0, yej,t>0; TopfV\y,0) = e„ yeVf. J 

In Eq. f l5.4ip ej is the standard Cartesian basis vector. 
Therefore 

B,ip,,t) = {B{{y,t))Y,, (5.42) 

b) If To = and ;Ui > then the solution of the stationary microscopic 
equations (15. 6p . (15.180 . and fl5.19p is given by formula 

V = ^ + Biiy).i-Wq + PfF), 

in which 

3 

Bi{y) = Y,U\y)^e,, 

i=l 

and the functions U^{y) are defined from the periodic boundary value prob- 
lem 

-piAU' + VR' = Bi, divyU' = 0, yeYf,] ^ ^ 

U' = 0, ye-f. J 

Thus 

B^ipi) = (B^((?/))y.. (5.44) 
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Matrices -Bi(/ii, t) and i?2(/^i) are symmetric and positively defined [3l Chap. 
8]. 

c) If To > and /ii = tlien in tlie process of solving the system fl5.6p . 
f l5.20p . and fl5.2ip we firstly find the pressure R{x,t,y) by virtue of solving 
the Neumann problem for Laplace's equation in Yf. If 

h{x, t) = -ToPf-^{x, t) - V^q{x, t) + pfF{x, t), 

then 

3 

R{x, t,y) = ^ Ri{y)ei (g) h{x, t), 

i=l 

where R^{y) is the solution of the problem 

AyRi = 0, ye Yf, VyRi n = n ■ Ci, G 7. (5.45) 

Formula f l5.39p appears as the result of homogenization of Eqs. f l5.20p and 

3 

S3 = 5^(Vi?.(2/))y,®e,, (5.46) 

i=l 

where the matrix (ml — ^3) is symmetric and positively definite P, Chap. 
8]. □ 

§6. Proof of Theorem [211 

6.1. Weak and two-scale limits of sequences of displacement and 
pressures. 

I. Let /ii < 00 and one of the conditions (12. 4p or 02.51) holds true. Then 
on the strength of Theorems 12.11 and 13 . 1 1 we conclude that sequences {x^iu^}, 
{p^} and {g^} two-scale converge to x{y)W{x, t, y), P{x, t, y) and Q{x, t, y) 
and weakly converge in L'^{VLt) to , p and q respectively, and a sequence 
{u'^{x, t)}, where u'^{x, t) is an extension of iu^(a;, t) from the domain into 
domain f2, strongly converges in L'^iVLj') and weakly in L^((0,T); W^iyi)) to 
zero. 

II. If < 00 and conditions fl2.5p hold true, then due to estimates (12.31) 
and (12. 6p the sequence {a\u'^} converges strongly in Li^iytT) and weakly in 
L^((0, T); 11^2^(1])) to a function it, and the sequence {tt*^} converges weakly 
in L'^[VLt) to a function tt. 

III. If /ii = 00, Pi^^rj^^ < 00 and < Ai < cxd, then on the strength 
of part 2) of Theorem 12.11 the sequences {a^e'^x'^'"^^}) {P^}^ {^^} and {q^} 
two-scale converge to functions x(2/)^(^) ?/)) Pi^^tjU)} n(a;,t,7/) and 
Q{x,t,y) and weakly in L'^{Qt) to functions , p, it and q respectively. 
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and the sequence {a^e ^u^} strongly converge in L'^{Qx) and weakly in 
L\{0,T);W^{n)) to the function u. 

o 

As before in previous section §5, we conclude that u E L^(0,T; W2 {^))- 
6. 2. Homogenized equations. 

I. If /ii < cxD and one of the conditions (12 ■4p or (12.51) holds true, then, as 
in the proof of Theorem 12.21 we construct a closed system of equations for 
the velocity v = dw^ /dt in the liquid component and for the pressures p and 
q, consisting of the modifications of Darcy's law (I5.37I) - (I5.39I) and boundary 
condition (I5.40p . in which we have u{x, t) = 0, and of the equations 

dp X dp 

p + ^op^^-g^ = q, ^^^^^ " ^ ^ ^ > 0- (6-1) 

We entitle the above described systems as Problem Fi, F2, or F3 depend- 
ing on the forms of the matrices Bi, B2, or B^, having places in Darcy's 
laws. 

II. Let fj,i < 00 and condition (12.51) holds true. We observe that the 
limiting displacements in the rigid skeleton are equal to zero. In order to 
find a more accurate asymptotic of the solution of the original model, we use 
again the re-normalization. Namely, let 

Then new displacements satisfy the same problem as displacements before 
re-normalization, but with new parameters 

a^ce^^, ax —^1, a,- — > a^-a^^. 

Thus we arrive at the assumptions of Theorem 12.21 Namely, the limit- 
ing functions u{x,t), TT{x,t), Il{x,t,y), and U{x,t,y) satisfy the system of 
micro- and macroscopic equations (15. 5p . (15.161) . (15.101) and (I5.17p . in which 
the pressure q is given by virtue of one of Problems Fi- F3. The only dif- 
ference from already considered case is in micro- and macroscopic continuity 
equations, because this equation depends on the value r]2- These micro- and 
macroscopic continuity equations coincide with Eq. (15. 5p and Eq. (I5.10p if 
we put there rjQ = rj^. 

Hence for u{x, t) and 7r(a^, t) there hold true the homogenized momentum 
equation in the form 

= div^.{Ag : D(x, u) + 5^div^.w + B{q - (g + tt) • 1} + pF, xe Q, (6.2) 

continuity equation (15.230 . in which we have r^o = ri2, and the boundary 
condition (^JUij. 

The tensor Aq, the matrices Cq,Bq and B^ and the constants Oq are 
defined from Eqs. (I5.30p . (I5.32p - (l5.33p . in which we have rjo = 772 and Aq = 1. 
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III. If Hi = oo, Pi ,r]^ < oo and < Ai < oo then re-normalizing by 



we arrive at the assumptions of Theorem 12.21 when /xi = 1, tq = Aq = Ai. 
Namely, functions , p, vr and u satisfy the following initial-boundary value 
problem in Qt'- 



div^{AiAg : B{x, u) + 5^div^n + B^p - (p + vr) ■ 1} + pF = 0, ' 

dw-f du 
" 'dt 



B2{l)-{-Vp + pfF), 



dt 

1 1 f 

— p H vr + diVxW^ = (m — l)diva;n, 

^ Pi Vi 

— vr + Cq : B{x, u) + Oodiv^n + a\p = 0. 



> (6.3) 



As before, tensor Ag, matrices Cq, Bq and Bf and constants defined 
by formulas fl5.30p . (15.321) - f l5.34p . in which we have r/o = Vi = -^i- 

Note, that here z/q = 0. Therefore the state equation p + UQp^^dp/dt = q 
becomes p = q. 

The problem is endowed by the corresponding homogeneous initial and 
boundary conditions. 



□ 

§7. Proof of Theorem [231 



7.1. Weak and two-scale limits of sequences of displacement and 
pressures. 

On the strength of Theorem 12. 11 the sequences {p^}, {g^}, {tt^} and {lu^} 
are uniformly in e bounded in L'^{Qt)- Then there exist a subsequence from 
{e > 0} and functions p, vr, q, and w such that as e \ 

w'^w, p'-^p, q'-^q, vr^^vr ^^(I^t). (7.1) 

Moreover, since Ao,/io > then the bound (11.17P imply 

V^w^ — > V^w weakly in L'^{VLt)- (7.2) 

e\0 

Due to limiting relations (17. ip . (17. 2p and Ngutseng's theorem, there exist 
one more subsequence from {e > 0} and 1-periodic in y functions P{x, t, y), 
Il{x,t,y), Q{x,t,y), and W{x,t,y) such that the sequences {p'^}, {vr"^}, 
{g^} and {Vw^} two-scale converge as e \ respectively to P, 11, Q and 

V^W + VyW. 
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7.2. Micro- and macroscopic equations. 

In the present section we do not consider functions of time t, whicli re- 
normalize pressures. As we liave shown before, we can ignore all functions 
of time. 

Lemma 7.1. Two-scale limits of the sequences {p^}, {tt^}, {g^} and {Vw^} 
satisfy inYT = Yx (0, T) the following relations 

—U + (l-x)idiv^w + diVyW) = 0; (7.3) 
Vo 

1 u dP 

-P + Xidzv.,w + dzVyW) = 0, Q = P + ^ — - (7.4) 

ot 

diVy{xi^omx, ^) + ©(y, -^)) + (1 - x)Ao(D(x, w) + D(y, W))) (7.5) 

-Vj,(g + n) = 0. 

Lemma 7.2. The weak limits p, tt, q w satisfy in Qt the following system 
of macroscopic equations: 

— 71 + (1 — m)diVxW + l.diVyW)Y = 0; (7.6) 

—p + mdiVxW + {diVyW)Yf=^, q = P + (7.7) 
p* y / Ql 

d^w , dw dW 

^^P^ = d^v4fio{mB{x, —) + (D(y, ^))y,) + (7.8) 

Ao((l - m)D(a;, w) + (D(y, W))y,) - (g + 7r)l) + pF. 

Proofs of these statements are the same as in lemmas 15.11 - 15. 3[ 
7.3. Homogenized equations. 

Lemma 7.3. Weak limits p, n, q and w satisfy in Qt the following system 
of homogenized equations: 

roP^ + V{q + TT)-pF= (7.9) 
duo 

diVx (A2 : D(x, + A3 : D(x, w) + B^diVxW + 

/ [Ki{t — t) : ©(x, w{x, r)) + B^{t — T)diVxW{x, T))dT), 
Jo 

— p + mdiVxW= (7-10) 
P* 

— / {C2{t — t) : 3{x,w{x,t)) + a2{t — T)diVxW{x,T))dT, 
Jo 
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— Ti + (1 — m) div^w = (7-11) 
Vo 

— / [Cslt — t) : ©(x, w{x, r)) + a^lt — T)diVxW{x, T))dT, 
Jo 

i/ere A2, A3 and A4 - fourth-rank tensors, B4, B^^, C2 and C3 - matrices and 
02 and 03- scalars. The exact expressions for these objects are given below by 
formulas (^M)-UM>- 

Proof. Let 

dw 

Z{x,t) = fioB{x,—) - \oB{x,w), Zij = Si ■ {Z ■ Sj), z{x,t) = div^w. 

As usual we look for the solution of the system of microscopic equations 
(E3D-(E5D in the form 

W= [ [W°(i/,t-r)^(a;,r)+ 5^iy^^(?/,t-r)%(a;,r)]rfr, (7.13) 
Jo 



.t 3 



P = X 



[ [P%y,t- t)z{x,t) + V P'^{y,t- T)Z,,{x,T)]dT, (7.14) 



Q = xiQoiy) ■ z{x,t) + J2 Qoiy) ■ z.^i^^t) + (7-15) 

/ [Q\y,t- t)z{x,t) + V Q'^{y,t- r)%(a;, r)]rfr), 

U={l-x)[ [U\y,t-T)zix,r)+y2w^iy,t-r)Z,,ix,T)]dT, (7.16) 

where 1-periodic in y functions P°, P^\ Qo, Q°, Q'^ , Q^^ 

n", n*-' satisfy the following periodic initial-boundary value problems in the 
elementary cell Y: 
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Problem (I) 



div^(x(//oD(y,^^) + 

(1 - x)(AoB(?/, W'^) - ((1 - + XQ'')^)) = 0; (7.17) 

1 u r) P^-^ 

-P'^ + X^^^yW^ = 0, Q'' = p^i + A 
p* p* at 

Liiij + (1 _ x)diyyW = 0, W'^y, 0) = W'i{y)- (7.18) 

■'70 

div,(x(/xoD(l/, W-^') + J^^' - Qo'D) = 0, (7.19) 

xiQ'i + ^odiv.iy^^) = 0. (7.20) 



Problem (//) 



Then 



diVj,(x(//o»(z/,^^) + 

(1 - x)(AoD(y, W^°) - ((1 - X)n« + XQ')1)) = 0; (7.21) 

x(-P° + div,W° + l) = 0, QO = pO + ^^; (7.22) 

(1 - x)(-n° + div,W-° + 1) = 0; (7.23) 

W'{y,0)^W'M, div,(x(//oB(y,T^S)-QoI))=0, (7.24) 

x(Qo + i^o(div^W° + 1)) = 0. (7.25) 



3 

A2 = Horn ® + /xqA^, 

3 



A3 = Ao(l - m) J] J^^' - AqAJ + //oAf(O), 

A4(i) = /^o^A{(t) - AoAf(t); (7.27) 



A{W = J]((/^oO(|/,^^(2/,t)))y, + (AoD(2/,W'^(2/,t)))yJ® J'^; (7.28) 
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B,{t) = /xo(ro(l/, ^)y, + Ao(ro(i/, W%y, t)))yj (7.29) 

3 

C2{t) = -C,{t) = Y,{dwyW'^{y,t))Y,r'; (7.30) 

a^(t) = -a,{t) = (div,iy°(y,t))y^,, B, = fiomy,W',{y)))Y,. (7.31) 

□ 

Lemma 7.4. Tensors A2- A4, matrices B4, B^, C2 and C3 anc? scalars 02 
anc? 03 are well- defined and infinitely smooth in time. 

If a porous space is connected, then the symmetric tensor A2 is strictly 
positively defined. For the case of disconnected porous space (isolated pores) 
A2 = and the tensor A2 becomes strictly positively defined. 

All these objects are well-defined if Problem (/) and Problem(//) 
are well-posed. The solvability of above mentioned problems and smoothness 
with respect to time follow, due to linearity, from the standard a'priory 
estimates. 

The symmetry of A2 proves in the same way as the symmetry of Aq. If the 
porous space is disconnected, then the problem fl7.19p has a unique solution 
linear in y, such that 

xmy,W^) + J'^) = 0. (7.32) 

The last equality implies A2 = 0. 

In this case the tensor A3 becomes strictly positively defined. Indeed 

3 



On the other hand, coming back to f l7.17p at initial time moment we see that 

{xfioHy, -g^(y^ 0)) ■ ^(^' ^o))y = 

-XoixHy, win : D(y, w'o'))y - i-u^' ■ n'=Vk=o. 

^0 

Moreover, due to (17.320 

(x/ioD(2/,^(2/,0)) : ©(y, W^o'')) = -(xD(y, ^(y, 0)) : 
which proves our statement. 
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